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ABSTRACT 


An  experimental  investigation  was  conducted  to  examine 
convection  heat  transfer  in  a  duct  of  rectangular  cross-sec¬ 
tion,  having  both  a  straight  and  curved  flow  passage.  One 
wall  was  kept  at  a  constant  temperature,  with  the  opposite 
wall  being  adiabatic.  Heating  of  the  air  flowing  in  the 
channel  was  accomplished  through  Joulean  heating  of  a  wall 
composed  of  Temsheet. 

The  experiments  were  conducted  for  steady  state  in 
Laminar  and  Transition  flows.  Development  and  progagation  of 
the  Taylor-Gortler  vortices  was  shown  to  enhance  the  heat 
transfer  rate  in  the  curved  section  as  compared  to  the 
straight  section. 
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I.  INTRODUCTION 


A.  TAYLOR-GORTLER  VORTICES 

Since  the  early  part  of  this  century,  considerable  data 
and  numerous  studies  have  demonstrated  that  fully  developed 
laminar  flow  along  a  concave  wall  does  not  remain  two-dimen¬ 
sional  [Refs.  1,  2,  3],  Instead,  the  flow  forms  a  system  of 
spiral  vortices,  of  counter  rotating  pairs,  whose  axes  are 
aligned  in  the  direction  of  the  principle  fluid  flow.  This 
phenomenon,  known  as  Taylor-Gortler  vortices,  is  the  result 
of  the  variations  in  the  centifugal  forces  acting  on  the 
fluid  particles.  Schematics  of  this  type  of  motion  are  shown 
in  Figures  1  and  2. 

In  a  channel  that  is  curved  in  the  streamwise  direction, 
the  fluid  particles  near  the  center  of  the  flow  cross-section 
are  subjected  to  higher  centifugal  forces  than  those  slower 
moving  particles  nearer  to  the  boundary  wall.  Consequently, 
the  tendency  is  for  the  fluid  in  or  near  the  center  of  the 
channel  to  move  outwardly  towards  the  concave  wall.  This 
action  then  causes  the  fluid  particles  near  the  boundary  wall 
to  move  in  the  spanwise  direction,  and  finally  move  radially 
inward,  replacing  the  outwardly  moving  particles.  Once  in  the 
center  of  the  channel  these  fluid  particles  come  under  those 
same  stronger  forces  of  the  free  stream  velocity  and  continue 
the  rotational  effect.  This  cyclic  motion  causes  the 


Schematic  of  Taylor  vortices  between  cylinders 


formation,  and  propagation,  of  th.e  counter-rotating  Taylor- 
Gortler  vortices-. 

Although  the  Taylor-Gortler  vortices  are  primarily  a 
laminar  flow  phenomenon  they  may  affect  the  transition  from 
laminar  to  turbulent  flows ,  and  are  therefore  important  to 
understand  {Ref.  2].  The  additional  mixing  provided  by  the 
secondary  motion  has  been  thought  to  account  for  the  increase 
in  the  rate  of  heat  transfer  from  concave  curved  walls  versus 
straight  walls  [Ref.  4] . 

The  Taylor-Gortler  vortices  appear  to  have  similarities 
to  other  vortex  flow  patterns  such  as  the  longitudinal  vor¬ 
tex  rolls  developed  in  the  laminar  forced  convection  heating 
of  fluid  layers  between  parallel  plates  [Ref.  5].  It  has 
also  been  thought  that  the  cross-hatching,  observed  in 
reentry  vehicles,  can  be  explained  at  least  in  part,  by  the 
presence  of  streamwise  vortices  [Ref.  &]  . 

There  are  many  possible  applications  that  could  result 
from  a  more  thorough  understanding  of  the  Taylor-Gortler 
vortices  and  their  effect  on  heat  transfer  and  fluid  flow 
characteristics.  Applications,  such  as  improved  turbine 
blade  cooling  [Refs.  7,  8,  9],  heat  exchanger  designs,  as  well 
as  other  such  engineering  designs  and  capabilities  could  take 
advantage  of  these  vortices  and  improved  heat  transfer 
characteristics  of  the  flow  in  the  presence  of  these  vortices. 
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B.  HISTORY 

The  instability  of  an  inviscid  fluid  rotating  symmetri¬ 
cally  about  an  axis,  as  in  curved  flow,  was  first  considered 
by  Lord  Rayleigh  in  1916  [Ref.  10],  In  assuming  that  the 
fluid  was  non-viscous,  he  determined  that  a  necessary 
condition,  for  the  stability  of  an  inviscid  fluid  motion,  was 
that  the  circulation  increased  with  increasing  radius.  G.  I. 
Taylor,  in  1923,  [Refs.  1,  11]  expanded  the  earlier  works  of 
Rayleigh  by  an  extensive  analytical  and  experimental  study  of 
viscous  fluids.  His  investigations  focused  on  the  flow 
between  two  cylinders,  in  which  the  inner  cylinder  rotated 
while  the  outer  cylinder  remained  stationary.  Taylor 
ascertained  that  such  couette  flows  became  unstable  when  the 
value  of  the  dimensionless  Taylor  number,  defined  as: 


exceeded  a  critical  value  of  41.3.  In  the  above  equation, 
'd'  is  the  width  of  the  gap,  assumed  small  when  compared  to 
'R^',  the  radius  of  the  inner  cylinder,  ’U^'  is  the  periphal 
velocity  of  the  inner  cylinder,  and  'Re'  is  the  Reynolds 
number.  Taylor  determined  that  for  those  cases  in  which  the 
value  of  the  Taylor  number  became  greater  than  the  critical 
value,  a  secondary  motion  developes  and  the  Taylor  vortices 
form. 
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Instability  of  a  similiar  nature  is  also  observed  when  a 
viscous  fluid  flows  in  a  curved  channel  due  to  a  pressure 
gradient  acting  along  the  channel  wall.  This  problem  was 
considered  analytically  first  by  W.  R.  Dean  [Ref.  12J  in 
1928,  for  a  channel  formed  by  two  concentric  cylinders, 
where  the  radius  of  the  inner  cylinder  was  large  in  comparison 
to  the  small  spacing  between  the  inner  and  outer  cylinder 
walls.  Dean,  in  his  studies,  concluded  that  there  would  be 
an  initiation  of  the  flow  instability,  and  the  formation  of 
vortices,  (similiar  to  the  Taylor  vortices),  when  the  Dean 
number,  defined  as: 

De  =  Re 

achieved  a  value  greater  than  36.  In  the  above  equation,  the 
Reynolds  number  is  based  on  the  mean  velocity  of  the  unper¬ 
turbed  flow,  'd'  represents  the  channel  half-width,  and  ' R^ ' 
is  again  the  inner  cylinder  radius.  The  analytical  work  of 
Dean,  was  later  verified  by  W.  H.  Reid  [Ref.  13] ,  using  an 
approximate  numerical  solution. 

H.  Gortler  [Ref.  2],  in  1940,  studied  the  stability  of 
laminar  boundary  layer  profiles,  on  curved  walls  under  the 
influence  of  small  disturbances.  In  his  study,  Gortler  found 
that  these  disturbances  were  similiar  in  nature  to  those 
studied  by  G.  I.  Taylor,  which  had  led  to  the  fluid  flow 
instability  in  the  form  of  vortices.  By  approximate  numerical 
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calculations,  Gortler  concluded  that  only  on  the  concave 
boundary  walls  were  the  amplified  disturbances  produced,  and 
that  the  overall  flow  profile  appeared  to  remain  laminar  in 
nature.  In  1955,  Gortler's  approximate  results  were  verified, 
with  an  exact  solution,  by  G.  Hammerlin,  as  reported  by  H. 
Schlichting  [Ref.  14]  and  then  further  substantiated  in  a 
more  extensive  numerical  analysis  by  A.  M.  0.  Smith  [Ref.  3]. 
More  recently,  verification  of  these  numerical  solutions  have 
been  demonstrated  by  use  of  hot  wire  anemometry,  laser  doppler 
systems,  and  flow  visualization  techniques,  [Refs.  15,  16], 
and  in  1976,  Y.  Aihara  [Ref.  17]  conducted  a  non-linear 
analysis  of  the  Gortler  vortices. 

With  the  growing  interest  in  the  effects  of  the  secondary 
flows  associated  with  the  Taylor-Gortler  vortices,  numerous 
specific  studies  have  been  published  concerning  the  influences 
of  these  vortices  on  the  transfer  of  heat  in  both  the  laminar 
and  turbulent  flow  regimes.  F.  Kreith  [Ref.  4],  in  1955, 
studied  the  influence  of  heat  transfer  with  respect  to  the 
curvature  of  the  boundary  wall  for  fully  turbulent  flows, 
and  concluded  that  the  heat  transfer  from  the  heated  concave 
boundary  wall  was  considerably  higher  than  that  transferred 
from  the  convex  boundary  wall  of  the  same  curvature  and  under 
similiar  turbulent  flow  conditions. 

In  1965,  L.  Persen  [Ref.  18],  considering  the  special 
cases  of  verv  high  and  very  low  Prandtl  number  fluids, 
related  the  increase  in  heat  transfer  rate  from  a  curved  wall 


to  the  presence  of  the  Taylor-Gort ler  vortices.  His  con¬ 
clusions  were  that  the  effect,  of  the  development  of  the 
vortices,  was  to  increase  the  rate  of  heat  transfer  through 
the  boundary  layer. 

One  of  the  first  experimental  works  on  the  effects  of  the 
Taylor-Gortler  vortices  with  respect  to  enhanced  heat  transfer 
was  by  P.  McCormack,  et  al.,  [Ref.  19],  published  in  1970. 

With  the  data,  they  concluded  that  it  was  necessary  to  retain 
the  non-linear  terms  of  the  flow  equation  in  order  that  the 
complete  theoretical  solution  of  the  effects  of  the  Taylor- 
Gortler  vortices  on  heat  transfer  be  explained.  In  1977, 

R.  Kahawita  and  R.  Meroney  [Ref.  20],  furthered  McCormack's 
work  and  concluded  that  the  higher  order  terms,  as  well  as 
the  normal  velocity  components  of  the  main  flow,  became 
increasingly  more  significant  to  the  calculations  at  smaller 
wave  numbers. 

In  considering  the  forced  convective  heat  transfer  for 
fully  turbulent  flows,  between  horizonal  parallel  flat  plates, 
Y.  Mori  and  Y.  Uchida  [Ref.  5],  in  1966,  reported  the 
formation  of  longitudinal  vortex  rolls  with  an  axis  parallel 
to  the  flow  direction.  These  streamwise-axis  directional 
rolls,  similiar  to  the  Taylor-Gortler  vortices,  formed  when 
a  critical  value  of  temperature  difference  between  the 
plates  was  exceeded.  In  1971,  M.  Akiyama,  et  al.,  [Ref.  21], 
confirmed  the  formation  of  these  vortices  at  a  critical 
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Rayleigh  number  of  17Q8,  for  the  case  of  heating  from  the 
bottom  plate  surface. 

There  has  been  a  rather  limited  amount  of  published 
literature  regarding  the  flow  and  heat  transfer  in  curved 
channnels  of  rectangular  cross-section.  Much  of  the  pub¬ 
lished  literature  involves  the  development  of  numerical 
approximations  and  solutions  for  heat  and  mass  transfer  in 
curved  ducts  of  various  geometries  (most  notably,  that  of 
circular  and  triangular  ducts).  K.  Cheng  and  M.  Akiyama 
[Ref.  22]  developed  a  numerical  solution  for  forced 
convection  heat  transfer  with  laminar  flows  in  curved  channels 
of  rectangular  cross-section,  but  only  foi  small  aspect 
ratios,  (B/D;  where,  'B'  is  the  channel  width  and  'D'  is  the 
depth  of  fluid  flow) .  Using  the  fully  elliptical  forms  of 
the  transport  equations,  G.  Yee  and  J.  A.  C.  Humphrey 
[Ref.  23] ,  developed  a  numerical  solution  for  flows  and 
boundary  conditions  similiar  to  Cheng  and  Akiyama.  K.  Cheng, 
et  al.,  [Ref.  24],  in  1974,  using  numerical  methods  solved  the 
Graetz  problem  for  a  curved  channel  of  square  cross-section 
CB/D-1) . 

In  1963,  B.  S.  Petukhov  and  V.  N.  Popov  [Ref.  25],  obtained 
an  analytical  expression  for  the  Nusselt  number  for  the  flow 
of  an  incompressible  fluid  with  variable  properties  in  a 
circular  tube.  In  1976,  A.  A.  Shibani  and  M.  N.  Ozisik 
[Ref.  26]  using  matched  asymptotic  expansion  techniques  for 
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the  case  of  uniform  wall  temperature,  solved  the  heat 
transfer  problem  between  parallel  plates  with  turbulent  flow, 
for  a  wide  range  of  Prandtl  numbers. 

As  to  analytical  and  experimental  studies  of  the  parallel 
plate,  channel  flow  heat  transfer  problem,  Y.  Mori,  et  al . , 
[Ref.  27] ,  obtained  results  for  hydrodynamically  fully 
de\eloped  flows  with  constant  wall  heat  flux,  in  curved 
channels  of  square  cross-section.  W.  M.  Kays  and  E.  Y. 

Leung,  [Ref.  28],  experimentally  obtained  solutions  for 
turbulent  flow  heat  transfer  in  a  concentric  circular  tube 
annulus  with  fully  developed  velocity  profile  and  constant 
heat  rate  per  unit  length,  for  a  fluid  of  Prandtl  number  0.7. 
M.  Durao  [Ref.  29]  and  J.  Ballard  [Ref.  30] ,  reported  results 
for  large  aspect  ratio  (B/D  =  40)  rectangular  cross-sectional 
channels  with  laminar  flows.  P.  F.  Brinich  and  R.  W.  Graham 
[Ref.  31] ,  reported  results  for  turbulent  flows  in  a  rectang¬ 
ular  curved  channel  with  an  aspect  ratio  of  6,  for  the  inner 
wall  heated,  the  outer  wall  heated,  and  both  walls  heated. 

In  1979,  F.  D.  Haynes  and  G.  D.  Ashton  [Ref.  32],  in  studying 
heat  transfer  from  a  river  to  its  ice  cover,  reported  results 
for  rectangular  cross-section  channels  of  aspect  ratio  10, 
and  compared  them  to  the  data  obtained  by  Ashton  [Ref.  33] , 


IT.  INTENT  OF  THE  STUDY 


The  primary  purpose  of  this  study  was  to  investigate  the 
effect  of  the  Taylor- Gortler  vortices,  on  the  enhancement  of 
th.e  heat  transfer  rate  in  a  curved  channel  of  rectangular 
cross -section  and  large  aspect  ratio.  This  investigation 
began  with  flow  velocities  that  were  in  the  laminar  flow 
regime  and  continued  into  the  transition  region,  between 
laminar  and  turbulent  flow.  Further,  this  study  compared  the 
heat  ransfer  rate  obtained  in  a  curved  section  with  the  rate 
obtained  in  a  straight  section  of  identical  aspect  ratio. 

It  was  expected  that  the  presence  of  the  Taylor-Gortler 
vortices  would  enhance  the  heat  transfer  process  at  all 
velocities  in  the  flow  regimes  investigated.  The  cause  of 
this  would  be  the  secondary  flow  velocity  components  of  the 
vortices,  transporting  heated  fluid  from  the  concave  boundary 
wall  inward,  while  simultaneously  displacing  the  cooler  fluid 
particles  in  the  center,  causing  them  to  move  outwardly 
towards  the  heated  concave  wall. 

This  study  was  conducted  using  a  single  rectangular  cross- 
sectional  channel,  that  incorporated  both  a  straight  and  a 
curved  test  section.  The  results  obtained  in  the  straight 
section  were  compared  to  those  results  obtained  in  the  curved 
section,  at  the  same  flow  velocities,  in  an  effort  to 
determine  the  effects  of  the  Taylor-Gortler  vortices  on  the 
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transfer  of  heat.  Additionally,  the  results  of  both,  the 
straight  and  curyed  section  tests  were  compared  to  the 
experimental  results  of  Ballard  JRef.  30]  and  Durao  [Ref.  29]. 

Those  results  of  the  straight  section  that  fell  into  the 
laminar  flow  region,  were  also  compared  to  the  analytical 
solution  as  given  by  P.  McCuen,  et  al.,  [Ref.  35];  while 
those  in  the  transition  region  were  compared  to  the 
experimental  results  of  Kays  and  Lueng  [Ref.  28] ,  and  the 
analytical  results  of  the  Dittus -Boelter  equation  [Ref.  38] 
and  Eckert's  equation  [Ref.  39].  While  these  equations  are 
for  fully  turbulent  flows  in  circular  tubes  they  have  been 
used  to  compare  turbulent  flows  in  other  geometry  channels 
when  the  diameter  of  the  tube  is  replaced  by  the  hydraulic 
diameter  (wetted  permimeter  divided  by  the  cross-sectional 
area) . 

The  straight  section  results,  for  transition  flow,  were 
also  compared  to  the  analytical  solutions  of  Petukhov  and 
Popov  [Ref.  25]  and  Shibani  and  Ozisik  [Ref.  26],  the  former 
for  circular  tubes  while  the  latter  is  for  parallel  plates, 
both  with  fully  turbulent  flows. 

The  curved  section  results  for  flows  in  the  transition 
region  were  compared  to  the  experimental  results  of  Brinich 
and  Graham  [Ref.  31],  and  also  to  the  solutions  of  Petukhov 
and  Popov,  and  of  Shibani  and  Ozisik,  mentioned  above. 


III.  EXPERIMENTAL  WORK 


A.  DESCRIPTION  OF  THE  APPARATUS 

A  channel  of  rectangular  cross-section  built  as  described 
in  references  29  and  30,  was  used  to  achieve  the  objectives 
of  this  study.  The  channel  was  manufactured  from  two  0.635 
centimeter  thick  sheets  of  plexiglas,  separated  by  0.635 
centimeter  spacers  whose  inside  lengths  were  used  as  the 
sides  of  the  channel. 

The  channel,  as  shown  in  Figure  3,  was  composed  of  a 
straight  section,  122.0  centimeters  long,  followed  by  a 
curved  section  of  180  degree  arc.  The  interior  concave  wall 
of  the  curved  section  has  a  radius  of  curvature  of  30.5 
centimeters.  The  rectangular  channel,  with  an  aspect  ratio 
of  40,  was  0.635  centimeters  high  by  25.4  centimeters  wide, 
with  a  cross-sectional  area  of  16.13  square  centimeters  and 
a  wetted  perimeter  of  57.07  centimeters.  A  cross-sectional 
view  of  the  channel,  for  both  straight  and  curved  sections, 
is  shown  in  Figure  4. 

The  working  fluid  used  in  this  investigation  was  air,  at 
room  temperature,  which  entered  the  channel  via  an  entrance 
bell  constructed  of  Plexiglas,  which  was  connected  to  the 
straight  section.  The  entrance  bell  was  designed  and 
manufactured  in  accordance  with  ASME  nozzle  standards,  with 
an  elliptical  curved  base  on  a  major  axis  equal  to  ten  inches 
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and  a  minor  axis  of  one  inch.  In  an  attempt  to  prevent  the 
introduction  of  dust  and  other  foreign  matter  into  the  flow 
channel,  cheese-cloth  was  attached  to  the  opening  of  the 
entrance  bell  nozzle.  After  passing  through  the  channel, 
the  flow  was  directed  through  an  aluminum  exhaust  nozzle 
which  was  connected,  by  3/4  inch  PVC  piping,  to  a  Fisher  and 
Porter  Company  variable  area  flow  meter  model  10A3565A. 

The  rotameter,  with  a  Bead  Guide  Tube  number  FP-1-35-G-10 
and  a  1GSVGT  float  serial  number  68T60,  at  100%  full  scale 
had  a  flow  rate  of  0.793  cubic  meters  of  air  per  minute 
(28.0  standard  cubic  feet  of  air  per  minute  at  one  atmosphere 
and  70°F) . 

The  flow  of  air  was  drawn  through  the  channel,  piping  and 
rotameter  by  a  single  electrically  driven  Spencer  Turbo 
Compressor,  30  Hp.  at  3500  rpm,  rated  at  550  cubic  feet  per 
minute  at  70°F  and  one  atmosphere.  The  flow  rate  was  con¬ 
trolled  by  the  use  of  valving  both  at  the  compressor  and  at 
the  exit  of  the  rotameter. 

To  obtain  the  experimental  heat  transfer  data  for  this 
study,  two  test  sections  of  the  channel  were  constructed: 

(1)  a  straight  test  section,  29.2  centimeters  in  length, 
located  far  enough  downstream  to  ensure  hydrodynamically 
fully  developed  flow,  with  a  heated  test  section  area  of 
741.9  square  centimeters.  (2)  a  curved  test  section,  28.3 
centimeters  in  length,  subtending  an  arc  of  53.1  degrees, 
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located  in  the  lower  half  of  the  concave  curved  portion  of 
the  channel,  with  a  test  section  area  of  717.4  square 
centimeters . 

In  each  of  the  test  sections,  the  outer  wall  was 
modified  as  shown  in  Figures  5  and  6.  Temsheet,  a  carbon 
impregnated  porous  paper  with  a  property  of  uniform 
electrical  resistivity,  was  glued  to  the  interior  surface  of 
the  outer  wall  of  each  of  the  two  test  sections.  The  flow  of 
air  was  then  to  be  heated  by  Joulean  heating  of  the  Temsheet. 
As  the  electrical  resistance  of  Temsheet  is  not  constant,  but 
rather  slowly  variable  with  temperature,  a  precision  resistor 
with  an  electrical  resistance  of  2.078  ohms  was  connected  in 
series  with  the  Temsheet  to  allow  calculation  of  the 


instanteous  power  being  supplied. 

The  variables  to  be  measured,  in  this  investigation,  were 
(1)  the  inlet  flow  temperature  for  each  of  the  test  sections 
(Tj^) ,  (2)  the  outlet  flow  temperature  at  the  exit  of  each  of 
the  test  sections  (Tq^^)  >  (.3)  the  temperature  of  the  heated 
boundary  wall  of  each  of  the  test  sections  (Two)  »  (4)  the 
temperature  of  the  unheated  boundary  wall  of  each  of  the  test 


sections  (.T^^)  ,  (.5)  the  flow  temperature  at  the  exit  of  the 
channel  (TEXIT) ,  (6)  the  temperature  between  the  two  insula¬ 
tion  layers  for  each  of  the  two  test  sections  (Tj^g),  (7) 
the  voltage  across  the  precision  resistor  (Vp^) ,  (8)  the 
heater  voltage  (V^) ,  and  (9)  the  flow  rate  through  the 
rotameter  (Q) . 
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Insulation 

^Temsheet 


6.  Curved  Test  Section,  detailed  schematic 


The  temperature  measurements  required  ahoye,  as  well  as 
all  other  temperature  measurements  in  this  investigation, 
were  made  with  Copper-Constantan  thermocouples.  All  the 
thermocouples  had  been  previously  calibrated  using  a 
ROSEMOUNT  Communtating  Bridge  model  920A,  and  a  ROSEMOUNT 
Constant  Temperature  Bath  model  913A,  as  outlined  by  Ballard 
[Ref.  30]. 

A  total  of  sixty-four  Copper-Constantan,  glass  insulated, 
30  gauge  thermocouples  were  constructed  and  inserted  in 
specific  locations  to  facilitate  the  measurement  of  the 
required  temperatures.  All  the  thermocouples  were  wired  to  a 
Hewlett  Packard  Data  Acquisition  System,  model  2010c,  which 
enabled  an  automatic  printed  record  of  the  instantaneous 
reading  for  each  thermocouple  location.  All  the  thermocou¬ 
ples  were  also  connected  to  a  Kaye  Instrument  Ice  Point 
Reference,  model  K140-4,  serial  1780. 

In  each  of  the  test  sections,  five  of  the  thermocouples 
were  inserted  between  the  first  and  second  layers  of 
ARMAFLEX  22  insulation,  and  another  five  inserted  between  the 
second  and  third  layers.  Each  of  these  sets  of  five 
thermocouples  were  connected  in  parallel  to  facilitate  the 
recording  of  an  average  value  of  the  temperature  (TINS) 
between  the  respective  layers  of  insulation.  Additionally, 
sets  of  four  thermocouples,  also  connected  in  parallel  to 
indicate  average  temperature,  were  inserted  into  the  channel 
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at  the  following  three  locations;  l.i  the  entn  i:  ■'  of  the 
channel,  (2)  the  exit  of  the  siraigh :  .  ■■■-r  section,  and  (5) 
the  exit  of  the  curved  test  section.  At  ;  ach  or'  these  three 
locations,  the  thermocouples  were  s e.'  five  -entimeters 
apart  in  the  spanwise  direction  to  a  Ll--w  reading  of  the 
bulk  temperature  of  the  flowing  air. 

In  each  of  the  two  test  section-,  eight  thermocouples 
were  placed  in  direct  contact  with  the  Temsheet  through 
small  diameter  holes  drilled  into  the  Fiexiglas  outer  walls, 
with  the  thermocouple  beads  electrically  insulated  from  the 
Temsheet  with  ENMAR  Heat  Resisting  Glyceryl  Phthalate.  The 
final  sixteen  thermocouples  were  placed,  eight  to  a  test 
section  in  the  Plexiglas  unheated  boundary  wall,  opposite 
those  positioned  in  the  Temsheet.  A  detailed  sketch  of  the 
thermocouple  placement  is  provided  in  Figure  7.  The  average 
heated  outer  wall  temperature  (T  ) ,  was  provided  by  those 
thermocouples  in  contact  with  the  Temsheet.  The  average 
unheated  inner  wall  temperature  (Tw^j ,  was  provided  oy  tne 
tnermocoupies  embedded  in  the  Plexiglas  wall. 

Tne  thermal  insulation  consisisted  of  layers  of 
ARMSTRONG  ARMAFLEX  22  Sheet  Insulation,  a  flexible  foamed 
plastic  material,  on  the  heated  side  of  each  test  section, 
and  black  foam  insulation  one-half  inches  thick  elsewhere. 
Three  layers  of  the  ARMAFLEX  22  insulation  were  used,  each 
layer  had  a  thickness  of  0.636  centimeters,  and  covered  an 
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area  slightly  larger  than  the  respective  test  section. 

These  layers  of  insulation  were  used  in  the  compution  of 
heat  flux  losses  to  the  environment,  sample  calculations 
of  these  losses  are  given  in  Appendix  C.  The  rest  of  the 
channel  was  covered  with  the  one-half  inch  thick  Black 
foam  insulation  to  reduce  the  losses  to  the  surroundings 
through  the  plexiglas  walls  as  the  heated  air  transited 
the  channel.  The  insulation  was  held  in  place  by  the  use 
of  Adhesive  Heat  Resistant  ventilation  Duct  Tape  as  seen 
in  Figure  8.  The  entire  length  of  the  channel  as  well  as 
all  the  connections  between  the  channel  and  the  rotameter, 
were  sealed  with  General  Electric  Silicone  Rubber  Sealant 
Caulk  to  ensure  that  there  was  no  leakage  of  air  into  the 
channel  or  rotameter  downstream  of  the  entrance  bell 
nozzle. 

Two  aluminum  electrodes,  0.318  centimeters  thick,  were 
inserted  in  between  the  Temsheet  and  the  plexiglas,  in  each 
of  the  test  sections,  as  shown  in  Figure  4.  The  plates, 
with  the  Temsheet  and  aluminum  electrodes  were  assembled 
in  the  channel  with  plastic  teflon  screws  for  the  purpose 
of  electrical  insulation. 

B.  EXPERIMENTAL  PROCEDURES 

Experiments  were  conducted  with  the  volumeteric  flow 
rate  varied  incrementally  from  0.198  to  0.674  cubic  meters 


Test  Channel 


of  air  per  minute  to  correspond  to  particular  Reynolds 
numbers.  The  data  was  recorded  and  calculations  performed 
for  each  flow  rate.  The  experimental  procedures  followed 
were  the  same  for  both  the  straight  ana  curved  test 
sections.  The  straight  test  section  runs  served  as  a 
baseline  for  the  curved  test  section  and  allowed  tor 
comparison  of  the  results.  The  appropriate  electrical 
power  required  to  heat  the  Temsheet  boundary  wall  to 
approximately  50°L  was  determined  through  preliminary 
experiments.  This  temperature  was  choosen  to  allow  at 
least  a  twenty  degree  celisus  temperature  difference 
between  the  unheated  and  heated  boundary  wall  average 
temperatures  and  to  facilitate  achieving  steady  state 
conditions  within  a  reasonable  time  span.  Estimates  of 
time  to  achieve  steady  state  conditions  and  decisions  as 
to  data  collection  were  also  based  on  these  preliminary 
test  runs. 

The  determination  as  to  what  constituted  a  steady 
state  condition  required  careful  consideration  in  view 
of  the  multiple  input  variables  and  their  effects  on  the 
results.  The  final  criteria  for  steady  state  was  based  on 
three  input  variables:  (1)  flow  rate  percentage,  as 
recorded  from  the  rotameter;  (2)  heated  boundary  wall 


36 


temperature,  as  recorded  from  the  thermocouples  in  that 
area;  and  (3)  the  heater  voltage,  as  recorded  on  the 
voltmeter.  When  these  conditions  varied  by  less  that 
two  percent  over  a  ten  minute  period,  the  experiment 
was  considered  to  have  achieved  steady  state.  Under 
these  conditions  the  assumption  of  steady  flow  and 
constant  heat  flux  from  the  heated  wall  were  assumed 
valid. 

Originally,  the  time  required  to  reach  this  steady 
state  condition  was  found  to  be  in  the  range  of  two 
hours,  accordingly  all  runs  were  at  two  and  one-half 
hours  durations.  Data  was  taken  on  one  hour  intervals 
for  the  first  two  hours  and  then  at  ten  minute  intervals 
for  the  next  thirty  minutes.  Later  it  was  found  that  the 
higher  flow  percentages  for  the  curved  test  section, 
required  about  thirty  minutes  more  time  to  full  achieve 
steady  state,  so  test  runs  with  a  flow  percentage  of 
sixty  percent  or  greater  were  extended  with  three 
additional  ten  minute  interval  readings  for  a  total  time 
of  three  hours. 

The  instantaneous  values  for  all  the  temperatures 
were  automatically  printed  by  the  digital  recorder  in 
millivolts.  A  standard  thermocouple  table  was  used  in 
conjunction  with  previously  obtained  thermocouple 
calibration  data  to  convert  the  millivolts  to  degrees 


37 


Celsius.  The  volumetric  flow  rate  was  read  directly  from 
the  installed  rotameter  as  the  percent  of  full  scale  flow. 
The  instantaneous  power  supplied  was  calculated  using  the 
relationship : 

VH  VPR 

Qp  =  - 

rpr 

Where  the  precision  resistor  voltage  (VpR)  and  the  heater 
voltage  (V^)  were  recorded  in  volts  from  the  installed 
voltmeter  in  the  digital  Data  Acquistion  System,  and  the 
electrical  resistance  of  the  precision  resistor  (Rpp)  had 
been  previously  determined  to  be  2.078  ohms. 

As  the  area  of  the  curved  test  section  was  slightly 
less  than  the  area  of  the  straight  test  section,  (0.0718 
vs.  0.0742  square  meters),  the  heat  flux  from  the  heated 
wall,  of  the  respective  test  sections  was  different.  The 
room  temperature  entering  the  test  section  was  checked 
with  the  use  of  a  mercury  thermometer  located  in  the 
vicinity  of  the  apparatus. 

During  the  preliminary  test  runs  of  the  straight 
section,  the  thermocouple  array  recording  the  exiting 
temperature  of  the  flow  of  air,  was  set  at  three 
millimeters  from  the  heated  wall.  This  was  satisfactory 
at  the  lower  volumetric  flow  rates,  but  at  the  higher 
flow  rates  it  led  to  calculations  in  which  the  power 
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convected  to  the  air  was  greater,  by  approximately  fifteen 
to  twenty  percent,  than  the  power  being  supplied  to  the 
heater.  It  was  determined  that  the  thermocouples  were 
not  indicating  the  correct  bulk  temperature  of  the 
exiting  flow  of  air,  but  rather  a  higher  value  of  temp¬ 
erature,  this  then  resulted  in  an  erroneously  high  value 
of  convected  heat  transfer,  and  higher  Nusselt  numbers. 

To  correct  this,  the  system  was  allowed  to  run  until  the 
walls  of  the  downstream  curved  section  achieved  steady 
state.  Then,  calculating  the  losses,  by  use  of  the 
thermocouples  in  between  the  insulation  layers  of  the 
curved  test  section,  and  applying  these  losses  to  the 
entire  curved  section  area,  the  energy  loss  was  calculated. 
The  flow  exiting  the  channel  was  assumed  to  be  thoroughly 
mixed  and  therefore  the  temperature  as  recorded  by  the 
positioned  thermocouples,  at  the  exit  of  the  curved 
section,  was  used  as  an  accurate  representation  of  the 
average  bulk  exit  temperature.  With  the  final  exit 
temperature  and  the  losses  now  known,  the  temperature  of 
the  air  as  it  exited  the  straight  section  was  calculated. 

To  correct  this  condition,  the  thermocouples  were 
inserted  further  in  to  the  channel,  or  raised  closer  to 
the  outer  wall  until  the  millivolt  reading  was  compatible 
with  the  calculated  correct  temperature.  After  this 
procedure,  accomplished  at  a  flow  rate  of  0.595  cubic 
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centimeters  of  air  per  minute  (751  of  full  scale)  was 
completed,  additional  test  runs  were  conducted  at  0.555 
and  0.635  cubic  centimeters  of  air  per  minute  (70%  and 
80%  respectively) ,  without  adjusting  the  thermocouple 
positions,  to  check  the  accuracy.  In  each  case,  the 
temperature  was  within  0.2SoC  of  what  the  energy  balance 
indicated  it  should  be. 

It  was  later  discovered  that  the  curved  section  re¬ 
sults  had  the  same  problem  at  flow  rates  greater  then 
0.595  cubic  centimeters  of  air  per  minute.  Since  the 
channel  construction  is  such  that  it  does  not  allow 
further  downstream  data  collection,  and  since  the  power 
difference  was  less  than  four  percent  it  was  decided  to 
proceed  with  the  correlation  of  the  data  and  experimenta 
tion  and  note  the  discrepancies. 
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IV.  PRESENTATION  OF  DATA 


A.  ANALYSIS 

The  uniform  electrical  resistivity  of  the  Temsheet, 
approximated  a  constant  heat  flux  surface  for  the  heated 
wall.  With  insulation  on  the  unheated  wall  the  heat  losses 
through  that  wall  were  negligible  and  the  wall  could  be 
considered  adiabatic.  With  the  large  aspect  ratio  of  40, 
the  channel  was  considered  to  be  approximately  parallel 
plates.  The  straight  portion  of  the  channel  upstream  of 
the  straight  test  section,  was  of  sufficient  length  to 
ensure  that  the  flow  was  hydrodynamical ly  developed  for 
the  flow  velocities  investigated.  Based  on  these  assump¬ 
tions,  the  experimental  configuration  was  modeled  as 
forced  convection  between  parallel  plates,  for  hydro- 
dynamically  developed  and  thermally  developing  flows, 
subjected  to  a  constant  heat  flux.  The  boundary  conditions 
are  one  wall  at  a  constant  specified  heat  flux,  and  the 
other  wall  insulated. 

To  analyze  this  condition,  several  quanties  were 
defined  as  follows: 

The  average  heat  transfer  coefficient  between  the 
heated  wall  and  the  flow  of  air  was  defined  by  the  equation: 

Q  .  hApL  AT 
air 
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where,  '0  .  '  was  the  heat  convected  to  the  air,  'H-' 
was  the  average  heat  transfer  coefficient,  'Ap^'  was  the 
area  of  the  heated  Temsheet  for  the  section  being  tested, 
and  'AT'  was  the  difference  between  the  average  heated 
wall  temperature  (T  )  and  the  average  bulk  temperature 
of  the  fluid  (Tg) .  The  average  bulk  temperature  of  the 
fluid  was  defined  as  the  arithmetic  mean  of  the  fluid 
inlet  bulk  temperature  (T-^)  and  the  fluid  outlet  bulk 
temperature  (TqUT) .  The  delta  temperature  (aT)  could  then 
be  expressed  as: 

AT  =  T-  T  =  T  -  Hi _ 

wo  B  wo  2 

The  actual  heat  convected  to  the  air  (Q  .  >  was 

air ) 

calculated  separately  using  th?  expression: 

Q  -  m  c  (T  -  T  ) 

air  ^  air  OUT  IN 

where  'cpapr'  was  the  specific  heat  of  the  air  at  constant 
pressure  and  'm'  was  the  mass  flow  rate  of  the  air.  The 
mass  flow  rate,  is  further  defined  as  the  product  of  the 
volumetric  flow  rate  of  the  air  and  the  local  flow  density 
of  the  air  as  it  transits  the  rotameter.  The  volumetric 
flow  rate  of  the  air  was  corrected  to  standard  conditions 


for  the  rotameter  and  the  local  flow  density  was  calculated 


assuming  "perfect  gas"  behavior  at  constant  pressure  and 
utilizing  the  temperature  of  the  air  as  it  exited  the 
channel  CT^xi T ) * 

The  average  Nusselt  number  was  then  calculated  from  the 
average  heat  transfer  coefficient,  with  the  expression: 

R  d 

mr  =  - 

Kair 

in  which  ’d’  is  the  height  of  the  channel  and  'K  .  ’is  the 

air 

thermal  conductivity  of  the  air  evaluated  at  the  average 
bulk  temperature. 

The  Reynolds  number  was  then  calculated  for  each  test 
run  as  follows: 


m  d 


where  again  'm'  and  ’ d '  are  the  mass  flow  rate  and  height 
of  the  channel  respectively,  'A  '  is  the  cross-sectional 
area  of  the  channel,  and  ’  u '  is  the  dynamic  viscosity 

all 

of  the  air  evaluated  at  the  exit  flow  temperature. 

For  the  experiment  runs  utilizing  the  curved  section, 


the  Dean  number,  defined  as: 


was  also  evaluated.  In  this  equation,  'd'  is  again  the 
height  of  the  channel,  and  * R ^  *  is  the  radius  of  curvature 
of  the  unheated  wall  surface  of  the  channel  (the  convex 
wall)  . 

The  heat  losses,  by  conduction,  through  the  heated 
plate  and  the  insulation,  by  conduction,  were  also  cal¬ 
culated  for  each  test  section  utilizing  the  expression: 

AT 

INS 

Q  =  - 

lo  (  x  /K  )  Apt 
INS  INS 


With  'ATj^g'  being  the  difference  in  temperature,  as  record¬ 
ed,  oetween  the  first  and  second  layers  of  AkMAILEa  22 
insulation,  and  that  recorded  between  the  second  and  third 
layers.  fXj^s’  was  the  thickness  of  each  layer  of  insula¬ 
tion,  *  ^ INS  *  was  t*ie  'thermal  conductivity  of  the 
insulation,  and  TApL’  was  again  the  area  of  the  heated 
surface . 


In  order  to  verify  that  Q„.  ,  as  calculated  bv  the 
temperature  rise  in  the  air,  was  actually  the  convection 
heat  transfer  from  the  Temsheet  heated  surface,  an 
estimate  of  the  radiation  heat  transfer  was  also  calculated 
for  each  test  run.  The  eight  thermocouples  in  the  lower 
Plexiglas  surface  measured  the  temperature  of  that  surface 
and  allowed  the  computation,  utilizing  the  Stefan-Boltzman 
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law  of  radiation,  of  the  radiated  heat  transfer  (Q  )  with 
the  expression: 


^r  = 


ff  (T 


wo 


Twi> 


R 


R 


where  'cr  *  is  the  Stefen-Boltzman  constant,  and  'R  *  is  the 
total  resistance  to  radiation  heat  transfer  between  the 
surfaces  as  computed  by: 


^PL  ewo  fewi 


In  the  above  equation,  ’ «wo ’  is  the  emissivity  of  the 

heated  Temsheet  wall  surface  and  ' e  .  '  is  the  emissivity  of 

wi 

the  unheated  adiabatic  wall  surface  (Plexiglas) , 

A  sketch  of  the  control  volume  and  the  energy  balance, 
with  sample  calculations  for  one  test  run,  is  given  in 
Appendix  B. 


B.  RESULTS 

The  data  obtained  from  the  experimental  test  runs  were 
evaluated  utilizing  the  expressions  described  in  the 
previous  section,  and  as  illustrated  in  Appendix  B.  The 
results  of  the  major  variables  are  shown  in  tabular  form, 
in  Table  I  for  the  straight  test  section,  and  Table  II  for 
the  curved  test  section.  A  plot  of  the  average  Nusselt 
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number  verses  Reynolds  number,  is  given  in  Figure  9  for 
comparison  of  the  two  test  sections.  The  uncertainty  bands 
have  been  indicated  for  the  plotted  experimental  data  and 
a  sample  calculation  for  the  complete  error  analysis  is 
given  in  Appendix  A. 

As  shown,  the  results  indicate  an  increase  in  the 
rate  of  heat  transfer  with  increasing  Reynolds  number  for 
both  the  straight  and  curved  test  sections.  Additionally, 
the  heat  transfer  rate  of  the  curved  test  section  was 
higher  than  that  of  the  straight  test  section  for  each 
Reynolds  number  investigated. 

In  investigating  the  results  obtained,  and  considering 
those  Reynolds  numbers  below  2000  (1025  based  on  this 
channel  configuration)  as  laminar,  at  least  squares 
straight  line  correlation  of  the  data  results  in: 

0.67 

Flu  =  0.06  5  Red 

for  laminar  flows,  and 

—  0.58 

Nu  =  0.117  Red 

for  transition  regime  flows,  in  the  curved  test  section. 
With  correlations  of 

__  0.25 

Nu  =  0.983  Re  ^ 
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for  laminar  flows,  and 


—  0.56 

Nu  =  0.110  Red 

for  the  transition  flows  in  the  straight  test  section. 

In  the  laminar  flow  regime,  the  increase  in  heat 

transfer  was  approximately  fifteen  percent,  while  in  the 

higher  transition  flow  region  the  increase  was  approximately 

thirty  percent.  This  compares  favorably  to  the  earlier 

studies  by  Kreith  [Ref.  4],  in  which  he  reported  that  the 

heat  transfer  rate  along  concave  walls  increased  from 

twenty-five  percent  to  sixty  percent  for  Reynolds  numbers, 

4  6 

based  on  hydraulic  diameter,  between  10  and  10  ,  and 
Ballard  [Ref.  30] ,  who  found  an  eleven  percent  increase 
in  heat  transfer  rate  for  laminar  flows. 
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Figure  9.  Straight  vs.  Curved  Section  results  for  the  Present  Study 


TABLE  I 


SUMMARY  OF  STRAIGHT  TEST  SECTION  RESULTS 


r 

\ 

Red 

Qair 

(W) 

IT 

(W/m2°C) 

AT 

(°C) 

>Hi 

i 

[f. 

784 

39.42 

21.28 

26.7 

5.08 

' 

941 

40.14 

22.56 

24.0 

5.38 

i 

1099 

43.26 

23.11 

25.2 

5.51 

1257 

45.13 

23.93 

25.4 

5.71 

•4 

1414 

51.80 

25.89 

27.0 

6.18 

‘1 

1571 

50.44 

26.59 

25.6 

6.34 

- 

1729 

53.50 

28.61 

25.2 

6.83 

1886 

57.01 

31.32 

24.5 

7.47 

► 

2042 

59.24 

31.99 

25.0 

7.63 

« 

t 

2203 

61.65 

33.50 

24.8 

7.99 

■ 

2361 

62.97 

35.30 

24.0 

8.42 

2517 

67.11 

35.44 

25.5 

8.45 

2675 

69.14 

36.39 

25.6 

8.68 
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TABLE  II 


SUMMARY  OF  CURVED  TEST  SECTION  RESULTS 


Red 

De. 

Qair 

(W) 

K 

(W/m2°C) 

AT 

(°C) 

Nu 

782 

114 

41.19 

23.63 

24.3 

5.64 

937 

137 

47.18 

26 . 57 

24.7 

6.34 

1092 

159 

50.88 

28.56 

24.8 

6.81 

1172 

171 

53.71 

30.49 

24.5 

7.27 

1251 

183 

57.16 

32.32 

24.6 

7.71 

i  1406 

205 

58.11 

33.01 

24.5 

7.88 

1570 

229 

62.63 

35.12 

24.9 

8.38 

1720 

251 

69.41 

37.23 

26.0 

8.88 

1883 

275 

70.90 

38.29 

25.8 

9.13 

2038 

297 

76.55 

41.66 

25.6 

9.94 

2196 

320 

80.36 

42.72 

26.2 

10.19 

2351 

343 

83.30 

45.98 

25.2 

10.97 

1  2505 

366 

88.82 

46.87 

26.4 

11,18 

2659 

388 

94.67 

49.88 

26.4 

11.90 

V.  DISCUSSION  AND  CONCLUSIONS 

With  the  measurement  of  the  unheated  wall  and  the  air  flow 
temperatures,  it  was  indicated  that  there  was  a  negligible 
difference  between  the  fluid  bulk  temperature  and  the  tempera¬ 
ture  of  the  unheated  wall.  Additionally,  the  energy  balance 
and  sample  calculations,  provided  in  Appendix  B,  demonstrate 
that  the  radiated  heat  transfer  from  the  unheated  wall  sur¬ 
face  is  minimal.  Thus  it  was  assumed  that  the  heat  transfer 
to  the  flowing  air  was  solely  by  convection  from  the  heated 
Temsheet . 

The  high  aspect  ratio  of  the  channel  as  discussed  pre¬ 
viously,  provided  the  original  basis  for  the  assumption  that 
the  experimental  apparatus,  as  configured,  was  one  of  infinite 
parallel  plates.  The  experimental  data,  for  those  Reynolds 
numbers  in  the  laminar  flow  region,  tend  to  substantiate  this 
assumption,  in  that  they  approach  the  theoretical  limit  for 
average  Nusselt  number  of  5.385,  for  parallel  plates  with  one 
wall  heated  at  a  constant  heat  flux  and  the  opposite  wall 
adiabatic  [Ref.  37].  For  a  Reynolds  number,  based  on  hydraulic 
diameter,  of  1835,  the  average  value  of  Nusselt  number  obtained 
was  5.382. 

With  the  formulation  of  the  problem  complete,  the 
experimental  data  obtained  can  now  be  compared  to  the  analytical 
solutions  and  experimental  results  of  the  same  or  similar 
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problem.  That  problem  being,  flow  between  infinite  parallel 
plates  with  one  wall  at  constant  heat  flux  and  the  opposite 
wall  adiabatic. 

A  comparison  of  the  experimental  results  obtained  in  this 
study  and  the  experimental  results  of  Durao  [Ref.  29]  and 
Ballard  [Ref.  30]  (who  investigated  the  problem  solely  for 
laminar  flows),  are  shown  in  Figure  10,  for  the  straight  sec¬ 
tion,  and  Figure  11,  for  the  curved  section.  In  both  these 
studies,  as  in  this  investigation,  the  channel  has  an  aspect 
ratio  of  40,  and  the  experimental  procedures  were  similar. 

It  can  be  seen  that  the  results  obtained  in  this  study,  in 
the  laminar  flow  regime,  appear  to  be  in  agreement  with  the 
data  obtained  previously  for  both  test  sections.  Additionally, 
numerical  correlations  of  the  present  data,  as  given  earlier, 
for  the  laminar  flows  compares  favorably  to  the  numerical 
correlations  of  both  Durao  and  Ballard,  for  Reynolds  numbers 
greater  than  600,  as  given  in  Appendix  C.  In  the  two  investi¬ 
gations,  mentioned  above,  comparisons  were  made  to  the  analy¬ 
tical  studies  of  McCuen,  et.  al.  [Ref.  35],  for  heat  transfer 
between  infinite  parallel  planes  with  constant  wall  temperatures 
and  heat  flux.  Ballard  also  compared  his  data  to  the 
Worse-Schimdt  solution  [Ref.  36],  for  heat  transfer  in  tubes 
and  annular  passages.  In  both  the  Durao  and  Ballard  studies, 
the  experimental  data  compared  relatively  well  to  the  predicted 
solutions  of  MCuen,  et  al.  A  curve  of  the  predicted  solution 
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of  the  McCuen,  et .  al.  study  is  shown  in  Figure  10.  The  re¬ 
strictions  of  that  study  were  as  follows: 

1.  the  velocity  profile  is  fully  established. 

2.  the  fluid  transport  properties  and  density  are  assumed 
constant . 

3.  axial  heat  conduction  is  negligible. 

4.  viscous  energy  dissipated  is  negligible. 

5.  conditions  are  invarient  with  time  (steady  state). 

All  these  restrictions  have  been  accounted  for  in  the 

design  of  the  experimental  configuration  and/or  verified  in 
the  experimental  processes.  A  dimensionless  axial  length 
coordinate  defined  as: 


was  calculated  for  'x'  equal  to  29.2  centimeters  (the  length 
of  the  straight  test  section).  In  this  equation,  ’Pr  1  is 
the  Prandtl  number  of  air,  taken  as  0.7,  ' '  is  the  hydraulic 
diameter  (1.239  centimeters),  and  'Re  '  is  the  Reynolds  number, 
based  on  the  hydraulic  diameter. 

Based  on  the  values  of  x*,  it  was  determined  that  the 
experiments,  for  laminar  flows,  were  in  the  thermal  entrance 
region.  It  was  therefore  decided  to  integrate  the  expression 
of  the  local  Nusselt  number  as  given  by  McCuen,  et .  al.  [Ref. 
35]  : 
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for  the  small  values  of  x,  in  order  to  obtain  the  following 
expression  for  the  average  Nusselt  number: 

3 


Nu  =  - 


4(1 


?  Q  •  -(2/3) 

—  ln  (1  "  076 709T  C1  +  qh1)x*  } 


^WO 


-)x! 


In  this  expression,  Q  .  was  the  heat  convected  to  the  fluid 

W  1 

from  the  unheated  wall,  and  Q  was  the  heat  convected  to  the 

'wo 

fluid  from  the  heated  wall.  In  accordance  with  the  previous 
assumptions,  Q  .  was  considered  equal  to  zero,  and  Q  was 
considered  equal  to  Q  .  ,  and  the  average  Nusselt  number  could 

a.  1  l 

then  be  calculated,  and  the  expression  reduces  to: 


Nu  «  In  (1  - 


oT6?'0"9  5' 


,-(2/3) 


) 


This  value  of  average  Nusselt  number  was  plotted  for  values 
corresponding  to  an  x*  less  than  0.042,  as  this  represents  the 
maximum  value  for  which  the  analytic  solution  for  Nusselt 
number  is  valid.  This  value  of  x*,  (0.042),  for  x  equal  to 
29.2  centimeters,  equates  to  a  Reynolds  number  value  of 
approximately  800  ,  based  on  hydraulic  diameter  (or  411  based 
on  the  channel  height) . 

The  graphical  representation  shows  that  the  experimental 
results  of  this  study  as  well  as  those  of  Ballard  and  Durao 
plotted  above  the  analytical  solution.  This  difference  can 
probably  be  explained,  at  least  in  part,  by  the  difference  in 
geometries  used  in  the  studies,  and  by  the  limitations  inherent 
in  any  experimental  work.  The  side  wall  effects  and  the  in¬ 
ability  to  totally  account  for  all  the  heat  transfer  processes 


Comparison  of  present  data  with  Durao,  Ballard  and  McCuen  for  Laminar 


Curved 


and/or  losses  could  account  for  the  differences  of  approximately 
nine  to  ten  percent. 

In  the  comparison  of  the  experimental  results  obtained  in 
this  study  for  those  flows  not  in  the  laminar  regime,  the 
comparison  and  correlation  was  made  to  those  analytical  solu¬ 
tions  and  experimental  results  for  completely  turbulent  flows. 
For  the  most  part,  the  analytical  solutions  are  valid  for 
Reynolds  number  values  greater  than  10  (equating  to  a  Reynolds 
number  of  5128  in  this  study) ,  and  tube  flow  (or  other  geo¬ 
metries  based  on  hydraulic  diameter) . 

A  graphical  comparison  of  the  results  obtained  in  the 
straight  test  section  runs  with  the  analytical  solutions  and 
other  experimental  results  is  made  in  Figure  12.  The  plot 
shows  the  curve  of  heat  transfer  in  a  straight  tube  for  con¬ 
stant  wall  temperature,  as  given  by  the  Dittus-Boelter  equations 
[Ref.  38]: 

Nu  =  0.023  Re  0,8  Pr  0,4 

and  by  the  empirical  relation  as  given  by  E.R.G  Eckert  [Ref. 


39]: 


Nu  = 


0.0396  Re  0,75  Pr 


(1  - 


Re 


— 7 

T7UZ) 


This  analytical  expressions  when  extrapolated  down  show  some 
agreement  with  the  values  obtained  at  flow  rates  between  forty 
and  seventy  percent  of  the  rotameter  full  scale.  Unfortunately, 
those  values  obtained  at  the  higher  flow  rates  showed  less  of 
a  correlation  than  those  mentioned  above.  These  results, 
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obtained  for  volumetric  flow  rates  higher  than  seventy  percent 
may  have  been  affected  by  the  fact  that  they  were  near  the  top 
end  of  the  rotameter  scale,  and  therefore,  may  inherently  be 
somehwat  less  accurate  than  the  mid-scale  readings.  The  ther¬ 
mocouple  problem  discussed  previously,  may  also  have  contri¬ 
buted  to  the  decrease  in  correlation. 

When  extrapolated  down,  the  experimental  results  of  Kays 
and  Leung  [Ref.  28],  for  heat  transfer  in  annular  passages, 
with  one  wall  heated  at  a  constant  heat  rate,  and  r*  equal  to 
1.0,  (where  r*  is  the  annulus  radius  ratio  r^/r  ),  equivalent 
to  parallel  plates,  appears  to  correlate  well  with  the  present 
results.  This  extrapolation,  was  again  utilizing  a  least 
squares  straight  line  correlation,  as  given  in  Appendix  C. 

Also  shown  are  the  analytical  solutions  of  Petuhkov  and  Popov 
[Ref.  25] ,  for  turbulent  heat  transfer  in  round  tubes  given 
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I 

I 

The  analytical  expression,  given  by  Shibani  and  Ozisik  [Ref. 

26],  for  turbulent  heat  transfer  between  parallel  plates: 

I 

Nu.  -  12  ♦  0.03  Re3  Prb 

where,  J 

a  -  0.3  8  -  - — 

(3.6  ♦  Pr) 

b  *  0.33  ♦  0.5  e'0,6  Pr 

valid  for  0.1<Pr  <10,  and  10^<Re  <10b,  is  also  plotted. 

While  no  conclusions  can  be  drawn  from  the  data  obtained 
because  it  is  not  in  the  range  of  Reynolds  numbers  of  these 
studies,  it  is  worthy  of  note,  that  the  present  data  as  extra¬ 
polated,  and  the  analytical  and  experimental  works  of  others 
can  be  graphically  shown  to  be  close.  It  should  also  be  noted 
that  the  investigations  of  Hsu  [Ref.  34],  Ashton  [Ref.  33,  and 
Haynes  and  Ashton  [Ref.  32]  all  plotted  experimental  results 
between  the  Petuhkov  and  Popov  expression  and  that  of  Shibani 
and  Ozisik,  as  shown  graphically  in  the  Haynes  and  Ashton 
study  above. 

The  curved  section  data  is  plotted  in  Figure  13,  and 
compared  to  the  experimental  work  of  Brinich  and  Graham  [Ref. 

31]  for  flow  and  mass  transfer  in  a  curved  channel.  The 
analytical  solution  of  Petukhov  and  Popov,  and  Shibani  and 
Ozisik  are  also  plotted  as  a  reference.  A  least  squares 
correlation  as  given  in  Appendix  C  was  used  to  extrapolate 
the  data  of  Brinich  and  Graham.  The  accuracy  of  the  data 
points  used  in  this  correlation  are  subject  to  errors,  in  that 


the  actual  values  were  not  given  in  their  study  and  had  to  be 
taken  from  a  plot  of  Stanton  number  versus  Reynolds  number  as 
given  in  the  reference. 

Irrespective  of  the  relatively  small  differences  in  the 
results  and  the  correlations,  the  trends  appear  to  be  conclu¬ 
sive.  The  heat  transfer  rates  can  be  seen  to  increase  with 
increasing  Reynolds  number,  and  also  that  the  heat  transfer 
rate  in  a  curved  channel  is  measurably  increased  over  that 
rate  in  a  curved  channel  is  measurably  increased  over  that 
rate  in  a  straight  channel.  This  increase  in  heat  transfer 
rate,  appears  to  begin  with  the  development  of  the  Taylor- 
Goertler  vortices  and  continues  into  the  transition  and  tur¬ 
bulent  flow  regions.  This  result  also  appears  to  have  been 
reached  by  K.  V.  Dement ' eva  and  I.  Z.  Aronov  [Ref.  41],  who 
in  1978  ,  investigating  hydraulic  resistance  and  heat  transfer 
in  curvilinear  channels  of  rectangular  cross-section, 
concluded  that  there  was  an  enhancement  of  heat  transfer 
intensity  in  curved  rectangular  corss - section  channels  as 
compared  to  straight  pipes  and  channels. 
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resent  Pata  for  Transition  Flows,  Straight  Secti 


Reynolds  Number 

Data  for  Transition  Flows,  Curved  Section 


VI.  RECOMMENDATIONS 


The  amount  of  experimental  work  remaining  in  the  study  of 
heat  transfer  in  curved  channels  is  significant.  Additional 
experiments  in  the  transition  region  and  in  the  fully  turbu¬ 
lent  region  of  fluid  flow,  should  be  conducted  to  refine  the 
correlations  for  average  Nusselt  number.  With  this  work,  it 
would  be  necessary  to  design  and  construct  a  larger  experi¬ 
mental  configuration,  capable  of  higher  volumetric  flow  rate 
and  greater  temperature  differences  between  the  heated  surface 
and  the  fluid  inlet  temperature.  Additionally,  care  should 
be  taken  to  ensure  that  the  temperature  of  the  exit  flow,  of 
each  test  section,  is  accurately  measured.  This  may  be 
accomplished  by  the  use  of  a  mixing  cup  arrangement  at  the 
exit  flow  locations.  In  a  channel  designed  and  constructed 
as  the  one  in  this  study,  care  must  be  taken  to  allow  for  the 
removal  of  the  mixing  cup  at  the  straight  section  exit  location, 
to  ensure  that  turbulence  is  not  induced  in  the  flow  by  the 
thermocouples  prior  to  the  fluid  entering  the  curved  section. 

Also  in  the  way  of  temperature  measurement,  a  thermocouple 
should  be  installed  at  the  entrance  to  the  rotameter,  or  pre¬ 
ferably  some  better  device  or  means  of  determining  the 
volumetric  flow  rate.  This  would  ensure  a  more  accurate 
measurement  of  the  fluid  temperature  for  an  evaluation  of  the 
fluid  properties  in  the  mass  flow  calculation. 


=  .0200 

The  uncertainty  in  the  mass  flow  rate  was  then  calculated  as; 

^  =  .0200 

m 

and  the  Reynolds  number  uncertainty  was 

- 

\  (.0200)  2  +  C.01987)  2  +  (.00223)  2  +  (.01976)  2=  .0345 
d 

Red  =  2042  +  71 

The  uncertainty  in  the  Nusselt  number  was  also  calculated 

in  an  identical  manner,  using  the  equation: 
tr-T  _  Ed 


and  the  uncertainty  as 

)2  +  (|p)2  +  (S^il)2 

Nu  »  E  Kair 

It  was  therefore,  necessary  to  first  calculate  the  uncertainty 
for  E.  Using  the  same  approach, 

11  a  .  At 


dh  _  JrdQair,  2  .  ,  PL, 2  .  ,dAT, 2 

r  ^ 

where  AT  =  T  -  T„ 
wo  B 


■  .0026  and  S-I  ■  .0048 
The  uncertainty  in  Qair  was  calculated  by: 


Qair  *  mCpair  CTqut  -  Tin) 
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and 


Qair 
where 


dQair  J, dm,  2  .  ,dCPair^2  .  ,d(;TOUT  TIN),2 

— * —  ~  V  +  It - >  +  It - rr - > 

'  '  LPair  OUT  IN 


m 


and 


—  =  0.0200 

m 

d(T 


dC 


Pair 


"Pair 


=  .0042, 


OUT  '  TIN  ^ 


TOUT  'tin 


=  .0623 


These  calculations  result  in  the  following  uncertainties 
gfil  =  .0656 

—  =  .0658 


h 

dNu 

Nu 


=  .0688 


It  is  therefore,  the  uncertainty  in  the  difference  between 
outlet  and  inlet  temperatures  of  the  flowing  air,  that  is 
ultimately  the  major  source  of  average  Nusselt  number  uncer¬ 
tainty.  This  then  is  the  basis  for  the  recommendation  to 
better  ascertain  the  exit  temperatures  for  each  test  section. 
The  calculations,  as  shown  above,  were  for  the  651  straight 
section  run.  Similar  calculations  for  the  curved  section  at 
65%  yield  the  following; 
dRe 


—  =  .0345 


U 

Red  =  2038  +  70 


dNu 

Nu 


.0563 
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Nu  =  9.34  +  0.525. 

The  values  for  the  uncertainties  for  all  the  variables 
the  curved  section  test  run  are  as  follows: 


QUANTITY 

UNCERTAINTY 

apl 

.0026 

Ac 

.0198 

CP 

.0042 

d 

.0199 

De 

.0398 

FT 

.0526 

Kair 

.0003 

Nu 

.0563 

• 

m 

.0200 

Q 

.0200 

Qair 

.0523 

Re 

.  0345 

tb 

.0073 

tin 

.0029 

tout 

.0106 

T 

W 

.0060 

T  -T 

OUT  IN 

.0482 

AT 

.0047 

U 

.0022 

P 

.0007 
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APPENDIX  B:  SAMPLE  CALCULATIONS 


A  sketch  of  the  control  volume  for  the  energy  balance, 
for  the  straight  section,  indicating  the  major  heat  transfer 
components  is  shown  in  Figure  14.  The  sample  calculations 
that  follow,  show  the  methods  and  equations  that  were  used 
to  obtain  these  heat  transfer  components.  The  heat  transfer 
components  were  required  to  satisfy  the  energy  balance,  and 
to  obtain  an  average  Nusselt  number  as  a  function  Reynolds 
number  of  the  flowing  air. 

The  curved  section  energy  balance  and  calculations  are 
similar 


FIGURE  14.  Energy  balance  in  straight  section 
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A.  SAMPLE  CALCULATIONS  DATA 


troom 

= 

21 . 6°C 

Q 

= 

55% 

VPR 

= 

3.2165  V 

VH 

= 

37.1210  V 

rpr 

= 

2.078 

T 

IN 

= 

0.839  mV  = 

21 . 82°  C 

TINS1 

= 

1.324  mV  = 

33 . 34° C 

TINS2 

1.583  mV  = 

39 .45°C 

Twol 

= 

1.840  mV 

Two  2 

= 

1.840  mV 

Two3 

= 

1.987  mV 

Two4 

= 

2.001  mV 

Two  = 

Two  5 

= 

2.112  mV 

Two6 

= 

2.125  mV 

Two  7 

3 

2.245  mV 

Two8 

= 

2.241  mV 

Twil 

= 

0.964  mV 

Twi2 

0.967  mV 

Twi3 

= 

1.012  mV 

Twi4 

3 

1.008  mV 

Twi  = 

Twi5 

3 

1.056  mV 

Twi6 

3 

1.047  mV 

Twi7 

3 

1.101  mV 

Twi8 

* 

1.100  mV 

tout 

3 

1.124  mV 

=  27 . 97°C 

texit 

— 

1.076  mV 

=  26 . 91°C 

apl 

= 

0.0742  m2 

50 .09°C 


25 . 76°C 


=  80.44° 
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A 

c 

= 

0.0016  m 

d 

= 

0.00635  m 

Kair 

= 

0.02662  w/m°C 

r 

^Pair 

uair 

= 

1.0057  KJ/Kg° C 
1.983  x  10'5  Kg/m 

kins 

= 

4.18  x  10  W/m° C 

axins 

s 

0.00635  m 

£ 

wo 

= 

0.70 

c 

"Wl 

3 

0.40 

F 

wo-wi 

a 

: 

1.0 

5.669  x  10"8  W/m2 

TEMPERATURE  CALCULATIONS 

1.  Temperature  Difference  in  Insulation  (AT^g) 

ATINS  =  TINS2  ’  TINS1  *  39,45  -  33.34  =  6.11°C 

2.  Bulk  Temperature  (Tg) 

Tb  ,  TIN  *  ..TOUI  ,  21  -82  *  27.97  .  24_90.c 

3 .  Mean  Temperature  Difference  (AT) 

AT  =  T  -  T-  -  50.09  -  24.90  =  25.19°C 
wo  B 


POWER  CALCULATIONS 


1. 


2. 


Power  Supplied  (Q^) 

«  .  VPR  VH  _  (3.2165) (37.121) 

QP  — C 2TUTF) - 


57.46  W 


Heat  Lost  Through  Outer  Plate  (Q^  ) 

q  -  ATINS  =  _ 6.11 _ _ 

Lo  AXINS/KINS'APL  (0. 00635) /( (4.1 8x10 _Z) (0.0742)) 


=■  2.98  W 
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ft 


3. 


Heat  Radiated  (Q  ) 
a.  Radiation  Resistance  (RR) 


1  -e 


n  =  WO 

R  APL£wo 


,  1  -e  . 

+  1  +  wi 

+  A  F  T  +  A  £ 
APLwo-wi  APL  wi 


+  _i_  -  1]  =  2^929 
APL  ewo  £wi  APL 


b.  Heat  Radiated  (Qr) 
a (T  4-T.  .4) 

WO  Wl  J 


Qr  = 


'K 


R, 


C5. 669x10 "8) (0.0742) (323. 254- 
2.929 


=  4.214  W 

Heat  Convected  to  Air  (Q  .  ) 

_ all 

a.  Corrected  Volumetric  Flow  Rate  (q) 

tfxtt(°r) 

q  =  q/0  =  q  *  rnrr 

Standard1 


(0.793)  (0.55) 
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540.13 

529.69 


b. 


=  0.0073mJ /sec 
Density  (p  ) 


298. 924) 


P 


101325 


c.  Mass  Flow  Rate  (m) 


m  =  qp  =  (0.0073) (1.177)  =  0 . 0086Kg/sec 
d.  Heat  Convected  to  Air  (Q  .  ) 

all 

Qair  =  m  CPair  (TOUT  '  TIN) 

=  (0.0086) (1.0057) (27.97-21.82) 

=  0.05319  KJ/sec 

=  53.19  W 


D.  AVERAGE  HEAT  TRANSFER  COEFFICIENT  (K) 

tr  _  Qair  _  53.19 

h  ■  aplAT  “  Co' on 2) (25.13") 

=  28.46  W/M2oC 


E.  AVERAGE  NUSSELT  NUMBER  (Nu) 


Nu 


Ed  _  (28.46) (0.00635)  _ 
Kair  (0. 02662) 


6.789 


F.  REYNOLDS  NUMBER  (Re) 

Re  =  m  d  =  (0.0086) (0,00635) 
d  yairAc  (1.983xl0'5) (0.0016) 


1721.19 


APPENDIX  C:  CORRELATIONS 


All  correlations  were  obtained  using  the  method  of  least 

squares  (first  degree  polynomial),  as  outlined  by  C.  F.  Gerald 

[Ref.  43] .  In  this  method  the  values  of  Reynolds  number  (x) 

and  Nusselt  number  (Y)  were  first  converted  to  the  appropriate 

natural  logrithemic  value,  after  which  the  summations  of  x. , 

2 

xi  ’  Yi>  anc*  xiYi  were  calculated.  The  values  of  these  quanti¬ 
ties  were  then  placed  in  a  matrix  to  solve  the  simultaneous 
equations : 

aoN  *  ajSXj  -  SYj 

aoIxi  *  ’  IxiYi 

where  N  is  the  number  of  data  points,  and  the  matrix  is 


N 

s*i’ 

a 

0 

’*Y1  ‘ 

2Xi 

al 

S 

5X1Y1 

After  attaining  the  values  of  aQ  and  a^,  the  equation  of 

the  line  became: 

Ln  Nu  *  a  +  a.LnRe 
o  1 

or 

_  a  a, 

Nu  *  e  Re 

The  quantities  for  the  various  studies  are  given  in  Table 


The  standard  deviation  for  the  present  study  curved  section 


correlation  was  calculated  as: 


where:  is  the  actual  value  of  average  Nusselt  number  ob¬ 

tained  in  the  test  runs 


y^  is  the  value  of  average  Nusselt  number  calculated 
from  the  correlation  equation 

N  is  the  number  of  data  points  used  in  the  correlation 
n  is  the  degree  polynomial  used. 

The  result  of  the  above  method  was 
a2  =  .02426 
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TABLE  III:  CORRELATIONS 


V 


! 


H 

O 

H 

O 

2 

HH 

2 

HH 

os 
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US 

E- 

US 
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< 
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< 

-J 

HH 

2 

HH 

2 

2 
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2 

^3 

so 

HH 

HH 

1— 1 

03 

2 

(2Q 

2 

o 

2 

2 

OS 
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< 

►j 
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< 

< 

2 

OS 

OS 

tu 
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j 

H 

E- 

H 

H 
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H 

LO 

LO 

>< 
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• 
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rH 

rH 
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• 
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to 

o 

LO 

a> 
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vO 

cn 
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'O 
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00 
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LO 
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-1 
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LO 

o 

00 

00 

rH 

X 

« 

• 

• 

• 

• 

• 

• 

w 

o 

to 

to 

C2 

vO 

to 

sO 

rH 

to 

rH 

CvJ 

o 

to 

V© 

OS 

to 

til 

E- 
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HH 

tO 

LO 

CM 

to 
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a> 
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o 

2 

eu 
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H 
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